Preclinical evaluation of the SARS-CoV-2 Mpro inhibitor RAY1216 shows improved pharmacokinetics compared with nirmatrelvir.

Although vaccines are available for SARS-CoV-2, antiviral drugs such as nirmatrelvir are still needed, particularly for individuals in whom vaccines are less effective, such as the immunocompromised, to prevent severe COVID-19. Here we report an α-ketoamide-based peptidomimetic inhibitor of the SARS-CoV-2 main protease (Mpro), designated RAY1216. Enzyme inhibition kinetic analysis shows that RAY1216 has an inhibition constant of 8.4 nM and suggests that it dissociates about 12 times slower from Mpro compared with nirmatrelvir. The crystal structure of the SARS-CoV-2 Mpro:RAY1216 complex shows that RAY1216 covalently binds to the catalytic Cys145 through the α-ketoamide group. In vitro and using human ACE2 transgenic mouse models, RAY1216 shows antiviral activities against SARS-CoV-2 variants comparable to those of nirmatrelvir. It also shows improved pharmacokinetics in mice and rats, suggesting that RAY1216 could be used without ritonavir, which is co-administered with nirmatrelvir. RAY1216 has been approved as a single-component drug named 'leritrelvir' for COVID-19 treatment in China.

Rational Design of Benzobisheterocycle Metallo-ß-Lactamase Inhibitors: A Tricyclic Scaffold Enhances Potency against Target Enzymes.

Antimicrobial resistance is a global public health threat. Metallo-ß-lactamases (MBLs) inactivate ß-lactam antibiotics, including carbapenems, are disseminating among Gram-negative bacteria, and lack clinically useful inhibitors. The evolving bisthiazolidine (BTZ) scaffold inhibits all three MBL subclasses (B1-B3). We report design, synthesis, and evaluation of BTZ analogues. Structure-activity relationships identified the BTZ thiol as essential, while carboxylate is replaceable, with its removal enhancing potency by facilitating hydrophobic interactions within the MBL active site. While the introduction of a flexible aromatic ring is neutral or detrimental for inhibition, a rigid (fused) ring generated nM benzobisheterocycle (BBH) inhibitors that potentiated carbapenems against MBL-producing strains. Crystallography of BBH:MBL complexes identified hydrophobic interactions as the basis of potency toward B1 MBLs. These data underscore BTZs as versatile, potent broad-spectrum MBL inhibitors (with activity extending to enzymes refractory to other inhibitors) and provide a rational approach to further improve the tricyclic BBH scaffold.

The Ableist stare: an interdisciplinary, narrative-driven exploration of staring at disabled bodies.

In this paper, I explore a phenomenon those with visible disabilities are all too familiar with: being stared at for their disabled bodies. Drawing on the interrelated fields of psychology, narrative, autoethnography and philosophy, I argue that staring at disabled bodies morally harms disabled people. This moral harm arises from the fact that not only does staring at disabled people fundamentally treat them as means to ends in which they cannot share, and thus, violates the Kantian formula of humanity, but also because this staring results in further, consequential harms for disabled people as well. In elaborating on these consequential harms, I draw largely on the works of disability ethicists Rosemarie Garland-Thomson and Elizabeth Barnes and argue that staring at disabled people contributes to the hermeneutical injustice disabled people face in their largely ableist world. Having identified these harms, I then explore the ameliorative potential of elevating disability narrative (with various disability narratives largely leading the discussion, including my own), drawing on Hilden Lindemann's Damaged Identities, Narrative Repair, and hope to call attention to the ways in which our broader structurally ableist world contributes to disabled people being stared at for their bodies in such harmful fashion.

Correction: 2-Mercaptomethyl-thiazolidines use conserved aromatic-S interactions to achieve broad-range inhibition of metallo-ß-lactamases.

[This corrects the article DOI: 10.1039/D0SC05172A.].

Ab initio quantum chemistry with neural-network wavefunctions.

Deep learning methods outperform human capabilities in pattern recognition and data processing problems and now have an increasingly important role in scientific discovery. A key application of machine learning in molecular science is to learn potential energy surfaces or force fields from ab initio solutions of the electronic Schrödinger equation using data sets obtained with density functional theory, coupled cluster or other quantum chemistry (QC) methods. In this Review, we discuss a complementary approach using machine learning to aid the direct solution of QC problems from first principles. Specifically, we focus on quantum Monte Carlo methods that use neural-network ansatzes to solve the electronic Schrödinger equation, in first and second quantization, computing ground and excited states and generalizing over multiple nuclear configurations. Although still at their infancy, these methods can already generate virtually exact solutions of the electronic Schrödinger equation for small systems and rival advanced conventional QC methods for systems with up to a few dozen electrons.

Enhancement by pyrazolones of colistin efficacy against mcr-1-expressing E. coli: an in silico and in vitro investigation.

Owing to the emergence of antibiotic resistance, the polymyxin colistin has been recently revived to treat acute, multidrug-resistant Gram-negative bacterial infections. Positively charged colistin binds to negatively charged lipids and damages the outer membrane of Gram-negative bacteria. However, the MCR-1 protein, encoded by the mobile colistin resistance (mcr) gene, is involved in bacterial colistin resistance by catalysing phosphoethanolamine (PEA) transfer onto lipid A, neutralising its negative charge, and thereby reducing its interaction with colistin. Our preliminary results showed that treatment with a reference pyrazolone compound significantly reduced colistin minimal inhibitory concentrations in Escherichia coli expressing mcr-1 mediated colistin resistance (Hanpaibool et al. in ACS Omega, 2023). A docking-MD combination was used in an ensemble-based docking approach to identify further pyrazolone compounds as candidate MCR-1 inhibitors. Docking simulations revealed that 13/28 of the pyrazolone compounds tested are predicted to have lower binding free energies than the reference compound. Four of these were chosen for in vitro testing, with the results demonstrating that all the compounds tested could lower colistin MICs in an E. coli strain carrying the mcr-1 gene. Docking of pyrazolones into the MCR-1 active site reveals residues that are implicated in ligand-protein interactions, particularly E246, T285, H395, H466, and H478, which are located in the MCR-1 active site and which participate in interactions with MCR-1 in ≥ 8/10 of the lowest energy complexes. This study establishes pyrazolone-induced colistin susceptibility in E. coli carrying the mcr-1 gene, providing a method for the development of novel treatments against colistin-resistant bacteria.

Bioisosteric Design Identifies Inhibitors of Mycobacterium tuberculosis DNA Gyrase ATPase Activity.

Mutations in DNA gyrase confer resistance to fluoroquinolones, second-line antibiotics for Mycobacterium tuberculosis infections. Identification of new agents that inhibit M. tuberculosis DNA gyrase ATPase activity is one strategy to overcome this. Here, bioisosteric designs using known inhibitors as templates were employed to define novel inhibitors of M. tuberculosis DNA gyrase ATPase activity. This yielded the modified compound R3-13 with improved drug-likeness compared to the template inhibitor that acted as a promising ATPase inhibitor against M. tuberculosis DNA gyrase. Utilization of compound R3-13 as a virtual screening template, supported by subsequent biological assays, identified seven further M. tuberculosis DNA gyrase ATPase inhibitors with IC50 values in the range of 0.42-3.59 µM. The most active compound 1 showed an IC50 value of 0.42 µM, 3-fold better than the comparator ATPase inhibitor novobiocin (1.27 µM). Compound 1 showed noncytotoxicity to Caco-2 cells at concentrations up to 76-fold higher than its IC50 value. Molecular dynamics simulations followed by decomposition energy calculations identified that compound 1 occupies the binding pocket utilized by the adenosine group of the ATP analogue AMPPNP in the M. tuberculosis DNA gyrase GyrB subunit. The most prominent contribution to the binding of compound 1 to M. tuberculosis GyrB subunit is made by residue Asp79, which forms two hydrogen bonds with the OH group of this compound and also participates in the binding of AMPPNP. Compound 1 represents a potential new scaffold for further exploration and optimization as a M. tuberculosis DNA gyrase ATPase inhibitor and candidate anti-tuberculosis agent.

Pyrazolones Potentiate Colistin Activity against MCR-1-Producing Resistant Bacteria: Computational and Microbiological Study.

The polymyxin colistin is a last line antibiotic for extensively resistant Gram-negative bacteria. Colistin binding to lipid A disrupts the Gram-negative outer membrane, but mobile colistin resistance (mcr) gene family members confer resistance by catalyzing phosphoethanolamine (PEA) transfer onto lipid A, neutralizing its negative charge to reduce colistin interactions. Multiple mcr isoforms have been identified in clinical and environmental isolates, with mcr-1 being the most widespread and mcr-3 being common in South and East Asia. Preliminary screening revealed that treatment with pyrazolones significantly reduced mcr-1, but not mcr-3, mediated colistin resistance. Molecular dynamics (MD) simulations of the catalytic domains of MCR-1 and a homology model of MCR-3, in different protonation states of active site residues H395/H380 and H478/H463, indicate that the MCR-1 active site has greater water accessibility than MCR-3, but that this is less influenced by changes in protonation. MD-optimized structures of MCR-1 and MCR-3 were used in virtual screening of 20 pyrazolone derivatives. Docking of these into the MCR-1/MCR-3 active sites identifies common residues likely to be involved in protein-ligand interactions, specifically the catalytic threonine (MCR-1 T285, MCR-3 T277) site of PEA addition, as well as differential interactions with adjacent amino acids. Minimal inhibitory concentration assays showed that the pyrazolone with the lowest predicted binding energy (ST3f) restores colistin susceptibility of mcr-1, but not mcr-3, expressing Escherichia coli. Thus, simulations indicate differences in the active site structure between MCR-1 and MCR-3 that may give rise to differences in pyrazolone binding and so relate to differential effects upon producer E. coli. This work identifies pyrazolones as able to restore colistin susceptibility of mcr-1-producing bacteria, laying the foundation for further investigations of their activity as phosphoethanolamine transferase inhibitors as well as of their differential activity toward mcr isoforms.

Interactions of hydrolyzed ß-lactams with the L1 metallo-ß-lactamase: Crystallography supports stereoselective binding of cephem/carbapenem products.

L1 is a dizinc subclass B3 metallo-ß-lactamase (MBL) that hydrolyzes most ß-lactam antibiotics and is a key resistance determinant in the Gram-negative pathogen Stenotrophomonas maltophilia, an important cause of nosocomial infections in immunocompromised patients. L1 is not usefully inhibited by MBL inhibitors in clinical trials, underlying the need for further studies on L1 structure and mechanism. We describe kinetic studies and crystal structures of L1 in complex with hydrolyzed ß-lactams from the penam (mecillinam), cephem (cefoxitin/cefmetazole), and carbapenem (tebipenem, doripenem, and panipenem) classes. Despite differences in their structures, all the ß-lactam-derived products hydrogen bond to Tyr33, Ser221, and Ser225 and are stabilized by interactions with a conserved hydrophobic pocket. The carbapenem products were modeled as Δ1-imines, with (2S)-stereochemistry. Their binding mode is determined by the presence of a 1ß-methyl substituent: the Zn-bridging hydroxide either interacts with the C-6 hydroxyethyl group (1ß-hydrogen-containing carbapenems) or is displaced by the C-6 carboxylate (1ß-methyl-containing carbapenems). Unexpectedly, the mecillinam product is a rearranged N-formyl amide rather than penicilloic acid, with the N-formyl oxygen interacting with the Zn-bridging hydroxide. NMR studies imply mecillinam rearrangement can occur nonenzymatically in solution. Cephem-derived imine products are bound with (3R)-stereochemistry and retain their 3' leaving groups, likely representing stable endpoints, rather than intermediates, in MBL-catalyzed hydrolysis. Our structures show preferential complex formation by carbapenem- and cephem-derived species protonated on the equivalent (ß) faces and so identify interactions that stabilize diverse hydrolyzed antibiotics. These results may be exploited in developing antibiotics, and ß-lactamase inhibitors, that form long-lasting complexes with dizinc MBLs.

Tautomer-Specific Deacylation and Ω-Loop Flexibility Explain the Carbapenem-Hydrolyzing Broad-Spectrum Activity of the KPC-2 ß-Lactamase.

KPC-2 (Klebsiella pneumoniae carbapenemase-2) is a globally disseminated serine-ß-lactamase (SBL) responsible for extensive ß-lactam antibiotic resistance in Gram-negative pathogens. SBLs inactivate ß-lactams via a mechanism involving a hydrolytically labile covalent acyl-enzyme intermediate. Carbapenems, the most potent ß-lactams, evade the activity of many SBLs by forming long-lived inhibitory acyl-enzymes; however, carbapenemases such as KPC-2 efficiently deacylate carbapenem acyl-enzymes. We present high-resolution (1.25-1.4 Å) crystal structures of KPC-2 acyl-enzymes with representative penicillins (ampicillin), cephalosporins (cefalothin), and carbapenems (imipenem, meropenem, and ertapenem) obtained utilizing an isosteric deacylation-deficient mutant (E166Q). The mobility of the Ω-loop (residues 165-170) negatively correlates with antibiotic turnover rates (kcat), highlighting the role of this region in positioning catalytic residues for efficient hydrolysis of different ß-lactams. Carbapenem-derived acyl-enzyme structures reveal the predominance of the Δ1-(2R) imine rather than the Δ2 enamine tautomer. Quantum mechanics/molecular mechanics molecular dynamics simulations of KPC-2:meropenem acyl-enzyme deacylation used an adaptive string method to differentiate the reactivity of the two isomers. These identify the Δ1-(2R) isomer as having a significantly (7 kcal/mol) higher barrier than the Δ2 tautomer for the (rate-determining) formation of the tetrahedral deacylation intermediate. Deacylation is therefore likely to proceed predominantly from the Δ2, rather than the Δ1-(2R) acyl-enzyme, facilitated by tautomer-specific differences in hydrogen-bonding networks involving the carbapenem C-3 carboxylate and the deacylating water and stabilization by protonated N-4, accumulating a negative charge on the Δ2 enamine-derived oxyanion. Taken together, our data show how the flexible Ω-loop helps confer broad-spectrum activity upon KPC-2, while carbapenemase activity stems from efficient deacylation of the Δ2-enamine acyl-enzyme tautomer.

Discovering Quantum Phase Transitions with Fermionic Neural Networks.

Deep neural networks have been very successful as highly accurate wave function Ansätze for variational Monte Carlo calculations of molecular ground states. We present an extension of one such Ansatz, FermiNet, to calculations of the ground states of periodic Hamiltonians, and study the homogeneous electron gas. FermiNet calculations of the ground-state energies of small electron gas systems are in excellent agreement with previous initiator full configuration interaction quantum Monte Carlo and diffusion Monte Carlo calculations. We investigate the spin-polarized homogeneous electron gas and demonstrate that the same neural network architecture is capable of accurately representing both the delocalized Fermi liquid state and the localized Wigner crystal state. The network converges on the translationally invariant ground state at high density and spontaneously breaks the symmetry to produce the crystalline ground state at low density, despite being given no a priori knowledge that a phase transition exists.

Postoperative Delirium Severity and Recovery Correlate With Electroencephalogram Spectral Features.

BACKGROUND: Delirium is an acute syndrome characterized by inattention, disorganized thinking, and an altered level of consciousness. A reliable biomarker for tracking delirium does not exist, but oscillations in the electroencephalogram (EEG) could address this need. We evaluated whether the frequencies of EEG oscillations are associated with delirium onset, severity, and recovery in the postoperative period. METHODS: Twenty-six adults enrolled in the Electroencephalography Guidance of Anesthesia to Alleviate Geriatric Syndromes (ENGAGES; ClinicalTrials.gov NCT02241655) study underwent major surgery requiring general anesthesia, and provided longitudinal postoperative EEG recordings for this prespecified substudy. The presence and severity of delirium were evaluated with the confusion assessment method (CAM) or the CAM-intensive care unit. EEG data obtained during awake eyes-open and eyes-closed states yielded relative power in the delta (1-4 Hz), theta (4-8 Hz), and alpha (8-13 Hz) bands. Discriminability for delirium presence was evaluated with c-statistics. To account for correlation among repeated measures within patients, mixed-effects models were generated to assess relationships between: (1) delirium severity and EEG relative power (ordinal), and (2) EEG relative power and time (linear). Slopes of ordinal and linear mixed-effects models are reported as the change in delirium severity score/change in EEG relative power, and the change in EEG relative power/time (days), respectively. Bonferroni correction was applied to confidence intervals (CIs) to account for multiple comparisons. RESULTS: Occipital alpha relative power during eyes-closed states offered moderate discriminability (c-statistic, 0.75; 98% CI, 0.58-0.87), varying inversely with delirium severity (slope, -0.67; 98% CI, -1.36 to -0.01; P = .01) and with severity of inattention (slope, -1.44; 98% CI, -2.30 to -0.58; P = .002). Occipital theta relative power during eyes-open states correlated directly with severity of delirium (slope, 1.28; 98% CI, 0.12-2.44; P = .007), inattention (slope, 2.00; 98% CI, 0.48-3.54; P = .01), and disorganized thinking (slope, 3.15; 98% CI, 0.66-5.65; P = .01). Corresponding frontal EEG measures recapitulated these relationships to varying degrees. Severity of altered level of consciousness correlated with frontal theta relative power during eyes-open states (slope, 11.52; 98% CI, 6.33-16.71; P < .001). Frontal theta relative power during eyes-open states correlated inversely with time (slope, -0.05; 98% CI, -0.12 to -0.04; P = .002). CONCLUSIONS: Presence, severity, and core features of postoperative delirium covary with spectral features of the EEG. The cost and accessibility of EEG facilitate the translation of these findings to future mechanistic and interventional trials.

ipie: A Python-Based Auxiliary-Field Quantum Monte Carlo Program with Flexibility and Efficiency on CPUs and GPUs.

We report the development of a python-based auxiliary-field quantum Monte Carlo (AFQMC) program, ipie, with preliminary timing benchmarks and new AFQMC results on the isomerization of [Cu2O2]2+. We demonstrate how implementations for both central and graphical processing units (CPUs and GPUs) are achieved in ipie. We show an interface of ipie with PySCF as well as a straightforward template for adding new estimators to ipie. Our timing benchmarks against other C++ codes, QMCPACK and Dice, suggest that ipie is faster or similarly performing for all chemical systems considered on both CPUs and GPUs. Our results on [Cu2O2]2+ using selected configuration interaction trials show that it is possible to converge the ph-AFQMC isomerization energy between bis(µ-oxo) and µ-η2:η2 peroxo configurations to the exact known results for small basis sets with 105-106 determinants. We also report the isomerization energy with a quadruple-zeta basis set with an estimated error less than a kcal/mol, which involved 52 electrons and 290 orbitals with 106 determinants in the trial wave function. These results highlight the utility of ph-AFQMC and ipie for systems with modest strong correlation and large-scale dynamic correlation.

Response to Comment on "Pushing the frontiers of density functionals by solving the fractional electron problem".

Gerasimov et al. claim that the ability of DM21 to respect fractional charge (FC) and fractional spin (FS) conditions outside of the training set has not been demonstrated in our paper. This is based on (i) asserting that the training set has a ~50% overlap with our bond-breaking benchmark (BBB) and (ii) questioning the validity and accuracy of our other generalization examples. We disagree with their analysis and believe that the points raised are either incorrect or not relevant to the main conclusions of the paper and to the assessment of general quality of DM21.

Virtual Screening Identifies Novel and Potent Inhibitors of Mycobacterium tuberculosis PknB with Antibacterial Activity.

Mycobacterium tuberculosis protein kinase B (PknB) is essential to mycobacterial growth and has received considerable attention as an attractive target for novel anti-tuberculosis drug development. Here, virtual screening, validated by biological assays, was applied to select candidate inhibitors of M. tuberculosis PknB from the Specs compound library (www.specs.net). Fifteen compounds were identified as hits and selected for in vitro biological assays, of which three indoles (2, AE-848/42799159; 4, AH-262/34335013; 10, AP-124/40904362) inhibited growth of M. tuberculosis H37Rv with minimal inhibitory concentrations of 6.2, 12.5, and 6.2 µg/mL, respectively. Two compounds, 2 and 10, inhibited M. tuberculosis PknB activity in vitro, with IC50 values of 14.4 and 12.1 µM, respectively, suggesting this to be the likely basis of their anti-tubercular activity. In contrast, compound 4 displayed anti-tuberculosis activity against M. tuberculosis H37Rv but showed no inhibition of PknB activity (IC50 > 128 µM). We hypothesize that hydrolysis of its ethyl ester to a carboxylate moiety generates an active species that inhibits other M. tuberculosis enzymes. Molecular dynamics simulations of modeled complexes of compounds 2, 4, and 10 bound to M. tuberculosis PknB indicated that compound 4 has a lower affinity for M. tuberculosis PknB than compounds 2 and 10, as evidenced by higher calculated binding free energies, consistent with experiment. Compounds 2 and 10 therefore represent candidate inhibitors of M. tuberculosis PknB that provide attractive starting templates for optimization as anti-tubercular agents.

QM/MM Simulations Reveal the Determinants of Carbapenemase Activity in Class A ß-Lactamases.

ß-lactam antibiotic resistance in Gram-negative bacteria, primarily caused by ß-lactamase enzymes that hydrolyze the ß-lactam ring, has become a serious clinical problem. Carbapenems were formerly considered "last resort" antibiotics because they escaped breakdown by most ß-lactamases, due to slow deacylation of the acyl-enzyme intermediate. However, an increasing number of Gram-negative bacteria now produce ß-lactamases with carbapenemase activity: these efficiently hydrolyze the carbapenem ß-lactam ring, severely limiting the treatment of some bacterial infections. Here, we use quantum mechanics/molecular mechanics (QM/MM) simulations of the deacylation reactions of acyl-enzyme complexes of eight ß-lactamases of class A (the most widely distributed ß-lactamase group) with the carbapenem meropenem to investigate differences between those inhibited by carbapenems (TEM-1, SHV-1, BlaC, and CTX-M-16) and those that hydrolyze them (SFC-1, KPC-2, NMC-A, and SME-1). QM/MM molecular dynamics simulations confirm the two enzyme groups to differ in the preferred acyl-enzyme orientation: carbapenem-inhibited enzymes favor hydrogen bonding of the carbapenem hydroxyethyl group to deacylating water (DW). QM/MM simulations of deacylation give activation free energies in good agreement with experimental hydrolysis rates, correctly distinguishing carbapenemases. For the carbapenem-inhibited enzymes, free energies for deacylation are significantly higher than for the carbapenemases, even when the hydroxyethyl group was restrained to prevent interaction with the DW. Analysis of these simulations, and additional simulations of mutant enzymes, shows how factors including the hydroxyethyl orientation, the active site volume, and architecture (conformations of Asn170 and Asn132; organization of the oxyanion hole; and the Cys69-Cys238 disulfide bond) collectively determine catalytic efficiency toward carbapenems.

In silico design of novel quinazoline-based compounds as potential Mycobacterium tuberculosis PknB inhibitors through 2D and 3D-QSAR, molecular dynamics simulations combined with pharmacokinetic predictions.

Serine/threonine protein kinase B (PknB) is essential to Mycobacterium tuberculosis (M. tuberculosis) cell division and metabolism and a potential anti-tuberculosis drug target. Here we apply Hologram Quantitative Structure Activity Relationship (HQSAR) and three-dimensional QSAR (Comparative Molecular Similarity Indices Analysis (CoMSIA)) methods to investigate structural requirements for PknB inhibition by a series of previously described quinazoline derivatives. PknB binding of quinazolines was evaluated by molecular dynamics (MD) simulations of the catalytic domain and binding energies calculated by Molecular Mechanics/Poisson Boltzmann Surface Area (MM-PBSA) and Molecular Mechanics/Generalized Born Surface Area (MM-GBSA) methods. Evaluation of a training set against experimental data showed both HQSAR and CoMSIA models to reliably predict quinazoline binding to PknB, and identified the quinazoline core and overall hydrophobicity as the major contributors to affinity. Calculated binding energies also agreed with experiment, and MD simulations identified hydrogen bonds to Glu93 and Val95, and hydrophobic interactions with Gly18, Phe19, Gly20, Val25, Thr99 and Met155, as crucial to PknB binding. Based on these results, additional quinazolines were designed and evaluated in silico, with HQSAR and CoMSIA models identifying sixteen compounds, with predicted PknB binding superior to the template, whose activity spectra and physicochemical, pharmacokinetic, and anti-M. tuberculosis properties were assessed. Compound, D060, bearing additional ortho- and meta-methyl groups on its R2 substituent, was superior to template regarding PknB inhibition and % caseum fraction unbound, and equivalent in other aspects, although predictions identified hepatotoxicity as a likely issue with the quinazoline series. These data provide a structural basis for rational design of quinazoline derivatives with more potent PknB inhibitory activity as candidate anti-tuberculosis agents.

Penicillanic Acid Sulfones Inactivate the Extended-Spectrum ß-Lactamase CTX-M-15 through Formation of a Serine-Lysine Cross-Link: an Alternative Mechanism of ß-Lactamase Inhibition.

ß-Lactamases hydrolyze ß-lactam antibiotics and are major determinants of antibiotic resistance in Gram-negative pathogens. Enmetazobactam (formerly AAI101) and tazobactam are penicillanic acid sulfone (PAS) ß-lactamase inhibitors that differ by an additional methyl group on the triazole ring of enmetazobactam, rendering it zwitterionic. In this study, ultrahigh-resolution X-ray crystal structures and mass spectrometry revealed the mechanism of PAS inhibition of CTX-M-15, an extended-spectrum ß-lactamase (ESBL) globally disseminated among Enterobacterales. CTX-M-15 crystals grown in the presence of enmetazobactam or tazobactam revealed loss of the Ser70 hydroxyl group and formation of a lysinoalanine cross-link between Lys73 and Ser70, two residues critical for catalysis. Moreover, the residue at position 70 undergoes epimerization, resulting in formation of a d-amino acid. Cocrystallization of enmetazobactam or tazobactam with CTX-M-15 with a Glu166Gln mutant revealed the same cross-link, indicating that this modification is not dependent on Glu166-catalyzed deacylation of the PAS-acylenzyme. A cocrystal structure of enmetazobactam with CTX-M-15 with a Lys73Ala mutation indicates that epimerization can occur without cross-link formation and positions the Ser70 Cß closer to Lys73, likely facilitating formation of the Ser70-Lys73 cross-link. A crystal structure of a tazobactam-derived imine intermediate covalently linked to Ser70, obtained after 30 min of exposure of CTX-M-15 crystals to tazobactam, supports formation of an initial acylenzyme by PAS inhibitors on reaction with CTX-M-15. These data rationalize earlier results showing CTX-M-15 deactivation by PAS inhibitors to involve loss of protein mass, and they identify a distinct mechanism of ß-lactamase inhibition by these agents. IMPORTANCE ß-Lactams are the most prescribed antibiotic class for treating bacterial diseases, but their continued efficacy is threatened by bacterial strains producing ß-lactamase enzymes that catalyze their inactivation. The CTX-M family of ESBLs are major contributors to ß-lactam resistance in Enterobacterales, preventing effective treatment with most penicillins and cephalosporins. Combining ß-lactams with ß-lactamase inhibitors (BLIs) is a validated route to overcome such resistance. Here, we describe how exposure to enmetazobactam and tazobactam, BLIs based on a penicillanic acid sulfone (PAS) scaffold, leads to a protein modification in CTX-M-15, resulting in irremediable inactivation of this most commonly encountered member of the CTX-M family. High-resolution X-ray crystal structures showed that PAS exposure induces formation of a cross-link between Ser70 and Lys73, two residues critical to ß-lactamase function. This previously undescribed mechanism of inhibition furthers our understanding of ß-lactamase inhibition by classical PAS inhibitors and provides a basis for further, rational inhibitor development.

Multiscale Simulations Identify Origins of Differential Carbapenem Hydrolysis by the OXA-48 ß-Lactamase.

OXA-48 ß-lactamases are frequently encountered in bacterial infections caused by carbapenem-resistant Gram-negative bacteria. Due to the importance of carbapenems in the treatment of healthcare-associated infections and the increasingly wide dissemination of OXA-48-like enzymes on plasmids, these ß-lactamases are of high clinical significance. Notably, OXA-48 hydrolyzes imipenem more efficiently than other commonly used carbapenems, such as meropenem. Here, we use extensive multiscale simulations of imipenem and meropenem hydrolysis by OXA-48 to dissect the dynamics and to explore differences in the reactivity of the possible conformational substates of the respective acylenzymes. Quantum mechanics/molecular mechanics (QM/MM) simulations of the deacylation reaction for both substrates demonstrate that deacylation is favored when the 6α-hydroxyethyl group is able to hydrogen bond to the water molecule responsible for deacylation but disfavored by the increasing hydration of either oxygen of the carboxylated Lys73 general base. Differences in free energy barriers calculated from the QM/MM simulations correlate well with the experimentally observed differences in hydrolytic efficiency between meropenem and imipenem. We conclude that the impaired breakdown of meropenem, compared to imipenem, which arises from a subtle change in the hydrogen bonding pattern between the deacylating water molecule and the antibiotic, is most likely induced by the meropenem 1ß-methyl group. In addition to increased insights into carbapenem breakdown by OXA ß-lactamases, which may aid in future efforts to design antibiotics or inhibitors, our approach exemplifies the combined use of atomistic simulations in determining the possible different enzyme-substrate substates and their influence on enzyme reaction kinetics.

Discovery of novel and potent InhA inhibitors by an in silico screening and pharmacokinetic prediction.

Aim: In silico screening approaches were performed to discover novel InhA inhibitors. Methods: Candidate InhA inhibitors were obtained from the combination of virtual screening and pharmacokinetic prediction. In addition, molecular mechanics Poisson-Boltzmann surface area, molecular mechanics Generalized Born surface area and WaterSwap methods were performed to investigate the binding interactions and binding energy of candidate compounds. Results: Four candidate compounds with suitable physicochemical, pharmacokinetic and antibacterial properties are proposed. The crucial interactions of the candidate compounds were H-bond, pi-pi and sigma-pi interactions observed in the InhA binding site. The binding affinity of these compounds was improved by hydrophobic interactions with hydrophobic side chains in the InhA pocket. Conclusion: The four newly identified InhA inhibitors reported in this study could serve as promising hit compounds against Mycobacterium tuberculosis and may be considered for further experimental studies.